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' A binomial-type operator on a stationary Gaussian process is in- 

troduced in order to model long memory in the spatial context. Con- 
' sistent estimators of model parameters are demonstrated. In particu- 

, lar, it is shown that In — d — Op{ ^^°^J^^ ), where d — (rii, d2) denotes 

• ' the long memory parameter. 

^ ' 1. Introduction and main results. A process obeying the spatial autore- 

gressive model 



(N 

> 



o 



(1.1) Xst = aXs-i,t + /3X.,t-i - a/3X,_i,t_i + e 



St 



was first studied by Martin [18], where — 1 < a,/3 < 1. Martin indicated that 
^S) . it is often desirable in practice that a process be reflection symmetric, that 

'nJ" is, pk£ = P-k-i = Pk-e = P-k,e for lags k and i, and that the autocorrela- 

tions have a simple form. These requirements led to the definition of model 
(1.1), which has autocorrelation p^^^ = qI'^I/JI^I for lags k and i. Spatial au- 
O ' toregressive models are shown by Tj0stheim [22] to be useful in studying 

geophysical quantities such as seismological data. Jain [14] indicates that 
these models can be applied to develop useful algorithms for image process- 
ing. Guiles and Gleeson [8] and Basu and Reinsel [2] illustrate the suitability 
of model (1.1) as an error term in a regression model used to analyze data 
collected in agricultural field trials. Moreover, empirical evidence of slow de- 
^ . cay of correlations between yield in two-dimensional agricultural field trials 

H I has received considerable attention (e.g., [10, 24, 25] and [20]). This led to 

the study of power law correlation functions by Whittle [25] and Besag [3] 
as an alternative to exponential decay. Martin [19] also indicates the impor- 
tance of long-range correlation structures in agricultural field experiments. 
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In addition, Professor Kwang-Yul Kim of the Department of Meteorology 
at Florida State University has communicated to us that many geophys- 
ical variables, such as ocean temperature, exhibit a well-extended spatial 
correlation structure (e.g., [16, 17]). The above references motivate the pos- 
sible need for inclusion of a long memory component when modeling certain 
spatial processes. 

It is assumed that the data-sites can be arranged in a square lattice. This 
commonly occurs in agricultural experiments where exactly one plant is 
located at site Tj0stheim [21] mentions that irregularly spaced data- 

sites can sometimes be replaced by a regular grid using the interpolation 
techniques of Delfiner and Delhomme [9] . 

Let C denote the field of complex numbers; define 

(p{zi,Z2,a,f3) = (1 - azi){l - j3z2), 

(1.2) 

i;{zi,Z2,a,(3) = [(l){zi,Z2,a,P)] ^ 

Suppose that BiXgt = Xs-i^t{B2Xst = Xs^t~i) is the backward shift operator 
on the first (second) index of Xgt'-, then (1.1) can be written compactly as 

(1.3) (l){Bi,B2,a,(3)Xst = est. 

It is assumed throughout this work that all processes considered are Gaus- 
sian. Asymptotic results on parameter estimators in the model (1.3) for both 
stationary and nonstationary cases can be found in [1, 4, 5, 15]. 

Our purpose here is to extend the work of Fox and Taqqu [11, 12] from 
time series to the spatial context by including a long memory component in 
the model (1.3). Memory in time series is modeled by use of the binomial 
operator (1 — B)'^, where d G i—^-, ^) denotes the memory parameter. The 
following operator is used in the spatial setting with two indices d = (^1,^2): 

(1.4) V'^ = (l-Bi)'^i 0(1-52)"^^ 

The operator V"' is defined by its corresponding power series representation 
in (21,22), that is, V^Xgt = Y^'k'£=oO'keXs-k,t-e, where the coefficients are 
found from the power series expansion of (1 — zi)'^^{l — Z'l)'^^ in the unit 
polydisc Ai(0) x Ai(0) with Ai(0) = {z G C: |z| < 1}. Let Z denote the set 
of all integers. Given ^ in (1.2), in (1.4) and the white noise process 
Est ~ WAA(0, o"^), the fractional autoregressive model of the form 

(1.5) 4>{Bx,B2,a,f5)SJ'^Xst = esu where s,t G Z and (fi,d2 G (-|, ^), 

is considered. It can be shown (see [6]) that the spectral density function of 
the stationary solution of (1.5) is 

^2 \l _ Q-ix\-2d-L _ ^-iyt-2d2 
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where 9 = {a,(],d), \a\ < 1, < 1 and \di\ < ^, i = 1,2. Moreover, the cor- 
responding autocovariance function when a = (3 = is 

^ (-l)^+^r(l-2di)r(l-2d,)a2 
' ' T{k-di + l)T{l-k-di)T{e-d2 + l)T{l-i-d2)' 

Let ^0 = (cto^Po^do) denote the true parameter value, where do = idio,d2o)- 
Recah that a second-order process {Xst :s,t £Z} is said to be stationary 
when E{Xst) = and, for each G Z, cov(X<j_|„fc , X^j) is independent of 
s,t gTj. Moreover, {Xgt -SjtGl^} is caUed strictly stationary provided all its 
finite-dimensional distributions remain invariant under translations. Denote 
Xn = 777 Em=i ^ke and define 

N 



E e'^^^+'y^-{XM-XN) 

kl=l 



to be the periodogram of the process. As suggested by Whittle [23] in the 
time series case, define ^at to be an argument 6 = {a,(3,d) making 

(1-7) ^%{e)=f i^^.dxdy 



ip f{x,y,0) 

a minimum, where /= [— 7r,7r], = [— 7r,7r] x [— 7r,7r] and f{x,y,9) is the 
spectral density function of the process. In our setting, it can be shown that 
a minimum exists when \a\ < r, \f3\ < r and \di\ < s, where < r < 1 and 
0<s<i, i = l,2. 

Some motivational comments concerning the estimator defined in (1.7) 
are in order. Assume that {Xst-s,t G Z} is a mean zero Gaussian pro- 
cess with spectral density function f{x,y,6) given in (1.6) and denote X = 
{Xii,Xi2, ... ,Xj\fN) ^ M{0,T]\f). Let Logx denote the natural logarithm of 
X. Then the log-likelihood function of X is 

Ln{9, a^) = - — Log 2^ - - Log jr^vl - 2^^^'^- 



Let Vat = (J-^Tn. Whittle [23] proves that \Vn{9)\ ^ I as N ^ 00 and, 
thus, LN{9,a'^) can be approximated by — ^Log27r(72 — 2^X'V^^(0)X 
when A'' is large. Fox and Taqqu ([11], page 518) use Parseval's identity to 
show that can be approximated by the more tractable Toeplitz matrix 
Bn = {hkimn). where = e^[('=-™)"+(^-")^]/-Hx, y, 0) dy, and this 

leads to the Whittle estimator given in (1.7). 

The primary results of this work are given below and proved in Sections 
2-4. The many details and statements whose proofs are not given here are 
all available on request from the authors. Our first result establishes strong 
consistency of the estimator Sjy. 



4 



BOISSY, BH ATTACH ARYYA, LI AND RICHARDSON 



Theorem 1.1. Assume that {X^jtSjtsZ} is a stationary Gaussian 
process having the spectral density function listed in (1.6), where \a\ < 1, 
\P\ < 1 and < di,d2 < ^. Then d]\f — > 9q almost surely as N ^ oo. 

The rate at which {9n} converges in probabihty is given in the next 
theorem. The notation -§gCrj^{Oo) is defined to be the partial derivative of 
o"^ evaluated at 6q. 

Theorem 1.2. Let {Xgt : s,t G Z} be a stationary Gaussian process hav- 
ing the spectral density function listed in (1.6), where \a\ < 1, \(3\ < 1 and 
< di , ^2 < ^ • Then: 

(a) N{9n - Oo) + NA]^^E{-§ga%{eo)) ^ J\f{0,U87r^^-^), where S and 
Aj\f are defined in (3.1), (3.2) and (3.5); 

(b) N{LogN)-^{9N-eo)^^^{6i3dw,624d2o,S33dio,6ud2o), whereK 
and 6ij are defined in (3.3). 

Observe that Aj\f in Theorem 1.2(a) depends on the abstract point 9n de- 
fined in (3.4). It fohows that Theorem 1.2(a) cannot be used for construction 
of a confidence set about Oq . It only shows that when 9^ — 9q is normalized 
by N, the second term in (a) provides just enough cancellation to produce 
a nondegenerate limit distribution. 

Remark 1.3. Many of the results given here are natural extensions of 
those given in the time series context by Fox and Taqqu [11]. However, a 
significant departure occurs with the bias term used in Theorem 1.2(a). Ac- 
cording to (4.3), the normalized convergence of the sequence X/|A;|<A'^(-^ — 
\k\/N)akiCki is required. The proof given by Fox and Taqqu ([11], page 529) 
shows that (in our notation) NfJ2\k\<Ni^ ~ \k\/N)akiCki — > and — > oo, 
for p <1. (They only used the case p = 1/2.) In the spatial case one has to 
deal with the asymptotic behavior of NE{-^a'j^{9o)), which amounts to con- 
sidering what happens to NJ2\k\<Ni^ ~ \t^\/^)(^kiCki as A^— > oo. Unlike in 
the time series setting [where the bias term VNJ2\k\<Ni^ ~ m/^)o-kiCki 
0] , this sequence goes to infinity. It is shown in this paper that (A/ Log'^ A^) x 
J2\k\<Ni^ — \k\/N)akiCki converges to a nonzero real number as A ^ oo. The 
point here is that the bias term, N Aj^^ E[-^a'j^{9o)], in the spatial context 
is of order Log^ A^ as A ^ oo. 

2. Proof of Theorem 1.1: Outline. A sequence of lemmas needed to es- 
tablish the strong consistency of estimators of model parameters is stated 
and the reader is referred to Boissy [6] for detailed proofs. The first lemma 
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gives a convenient method (in our setting) for finding the Fourier series 
representation of a spectral density function in terms of the orthogonal 
set ^e^(^^+^y) -.kjiGZ} with respect to the product Lebesgue measure A = 
Ai X A2 on I^. Let L'P{I'^) denote the set of all complex- valued functions for 
which \f\P d\ < 00, where p > 0. 

Lemma 2.1. Let g he a complex-valued function defined on the closed 
unit polydisc Ai(0) x Ai(0) in C x C for which (/(e"*^, e"^?') G L'^{P). Sup- 
pose that g is an analytic function defined on Ai(0) x Ai(0) with power series 
representation g{zi,Z2) = J2'k',£=o^kiZiZ2, and assume that X{{x,y) £ P:g 
is discontinuous at (21,2:2) = (e~*^,e~*^)} = 0. Then 5(e~'^,e~*^) has the 
Fourier series expansion J2'k',e=o^ki^~^^''^^^^^ ^^^2/ ^/Sfc^£=o l^'fc^P < 

00. 

The representation described in Lemma 2.1 can be used to verify Lemma 2.2(a). 

Lemma 2.2. Let f{x,y,0) denote the spectral density function listed in 
(1.6). Then: 

(a) feL^L^) and JpLog^f{x,y,9)dxdy = 0, 

(b) I12 fifll dxdy > 47r2 when 61^62. 

Lemma 2.3. Assume that {Xst-s,t£ Z} is a stationary process defined 
on the underlying probability space (Q,^,P) and having mean fi, autoco- 
variance 7 and spectral density function h{x,y,9Q). Suppose that Xj\f 
H and ■^J2^t=ii-^s+k,t+e — lj){Xst — A*) — > l{k,i) almost surely, and that 
(7 : (RP"'"^, !Bp+2) (R, 53) is bounded, Borel measurable and periodic with 
g{x + 27rfc ) = g{x) for each x G W^'^ and k G IP^'^ . If l'^ x K is a com- 
pact subset of which is contained in the set of all continuity points of 
g and In is the periodogram of the process, then Jj2 g{x,y,9)Ij\f{x,y) dxdy ^ 
Jj2g{x,y,6)h{x,y,6())dxdy uniformly in 9 £K, almost surely [P]. 

The omitted proof uses the technique of Hannan ([13], Lemma 1) in the 
time series case by uniformly approximating g{x,y,9) over P x K with an 
A^th order Cesaro sum (see [26], page 304, Theorem 1.20). Since a station- 
ary Gaussian process is strictly stationary, it can be shown that Xj\f — > ji 
and ^J2^^t=ii^s+k,t+e - lAi^st - IJ) ^ likj) almost surely and, thus. 
Lemma 2.3 is applicable in our context. 

Proof of Theorem 1.1. Combining Lemmas 2.2-2.3 above with the 
argument given by Hannan ([13], Theorem 1) for the time series case shows 
that 9n ^ 9q almost surely (see [6]). □ 
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Next, a sequence of estimators of that is more suited for computational 
purposes than (1.7) is given. Following Hannan ([13], page 133), define On 
to be the argument making 

(2.1) -alio)- ' ^ 



a minimum, where Wg = 2Trs/N. Under the hypothesis of Theorem 1.1 above. 



it can be shown that On — > almost surely. Once On has been found, d'j^{ON) 



in (2.1) can be used as a strongly consistent estimator of cj^. 

3. Proof of Theorem 1.2(a). In the time series context, Fox and Taqqu 
[12] used a combinatorial method in order to establish convergence to nor- 
mality of certain sequences of quadratic forms determined by Toeplitz ma- 
trices. Lemma 3.2 below allows one to extend the above result to the spatial 
context when the coefficient matrix can be expressed as a finite sum of Kro- 
necker products. This lemma is used to prove Theorem 1.2(a). 

Given heL^{P), aMmn = Si2e^^^^-'^>+^^-''^y^h{x,y)dxdy is a Fourier 
coefficient of h and TN{h) = (flfc^mn) denotes the corresponding N"^ x A^^ 
(block) Toeplitz matrix. Elements in {akimn) ^-re arranged lexicographi- 
cally beginning with row 1, followed by row 2, and so on. In particular, 
element akimn appears in row {k — 1)N + I and column {m — 1)N + n of 
TN{h). Recall that the Kronecker product of matrices A and B is defined by 
B = {ukeB), and the following properties are needed: nj=i(^i *^ ^j) = 
{Ylj=i Aj) (g) {llj=i Bj) and Tr(^ (g> B) = Tt{A) ■ Tr{B), when the matrices 
are compatible and Tr denotes the trace. 

Definition 3.1. A function /i:/— >]R is called admissible provided the 
following conditions are fulfilled: 

(A.l) h is symmetric and integrable; 

(A. 2) the set of discontinuities of h has Lebesgue measure zero; 

(A. 3) for each fixed 6 > 0, h is bounded on [5, vr]; 

(A. 4) there exists a < 1 such that h{x) = 0(|x|~") as x ^ 0. 

Given h,k:I the product of h and k is defined to be the function 

(h (Si k){x,y) = h{x) ■ k{y) for each x,y £ I. 

Lemma 3.2. Assume that fs = fsi fs2, Qt = On ^ gti, f = J2T=i fs, 
g = X]r=i 9t o-nd each fgj, gtj is admissible on I, j = 1,2. Moreover, suppose 
that fsj{x) = Odx]""^), gtj{x) = 0{\x\~^) as x ^ and a + P < ^. Then: 



(a) ^TT[TN{f)TN{9)? - (4vr2)3/^,[/g]2dxdy, 

(b) ^Tt[TN{f)TN{g)f^O whenk = ?,A..... 
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Proof. First, observe that if hst = fsi ® gt2, then T]\f{hst) = T/v(/si) ^ 
TN{9t2) = {akeTN{gt2))- Indeed, a typical element in {akiT^ {gt2)) is of the 
form aubran = Jj e'^''-'^> fsiix) dx-Jj e*(""-")f c/t2(?/) dy = e*[('=-^)^+("^-")j/] x 
hst{x, y) dx dy. The latter quantity appears in row N{k — 1) + m and column 
N{£ — 1) + n of both T]\f{fsi) Tj\f{gt2) and T]\f{hst), which establishes the 
equality. 

(a) Employing Theorem 1(a) of [12] and properties of the Kronecker prod- 
uct, 

^Tr[r^(/)r^(<7)]2 

= ]^ E T^[TN{fs)TN{gt)TN{fu)TN{gv)] 

s,t,u,v 

^ Tt[{TNifsl)'i^TN{fs2)){TN{gtl)^TN{gt2)) 



'',t,U,V 



X {TN{ful)^TN{fu2)) ® {TN{gvl)^TN{gv2))] 

N 



E 77 T^^[TN{fsl)TN{gtl)TN{ful)TN{gvl 



^ Yl (271")^ / [/si5ti/ni5t)i]'^a;- (27r)^ / [fs2gt2fu2gv2]dy 

= (47r2)3 f [fgfdxdy. 

The above application of Theorem 1 of [12] is valid since each fsj{gtj) has 
the same order as x ^ 0, respectively. It is not necessary that all /sj's be 
equal in the proof of Theorem 1 of [12]. 

(b) Suppose that k is an integer exceeding 2. First, assume that 
k(a + f3) < 1. An extension of the argument used in part (a) shows that 
^Tr[T,v(/)r;v(9)]'^ ^ (47r2)2'=-i/^,[/5]fcdxdy, and since k > 2, 
Tic[Tj\[{f)T]\f{g)]'' — > 0. The case when k{a + /?) > 1 is verified in a similar 
manner by employing Theorem 1(b) of [12]. □ 

Given the process {Xst-s,t G Z}, recall that = ■^J2^t=i-^st, and 

define X'j^ = {Xu — Xn,Xi2 — X^n, ■ ■ ■ ,Xi]^ — Xm,X2i — Xn, . . . ,X2n — 
Xat, . . . , Xjvi — Xjy, . . . , Xj\[j\f — Xj\f). Verification of the following result is 
based on the fact that the normal distribution is determined by its moments. 
The details of the proof are given in [6] . 
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Lemma 3.3. Let {Xgt ■ s,t £ Z} be a stationary Gaussian process hav- 
ing mean n and spectral density function f , and let g : (/^, *B2) — > (K, !B) be a 
bounded measurable function obeying g{—x, —y) = g{x, y) for each (x, y) G P . 
Moreover, assume that f = J2T=i fs [q = ELi 9t), fs = fsi ® fs2 (gt = Qti ® 
gt2) and each fsjiotj) is an admissible function with parameter a(/3) in 
Definition 3.1, where a -\- (3 < ^. If An = TN{g), then jjIX'j^AnXn — 
E{X'^AnXn)] ^ AA(0, 6^), where 6^ = UStt^ Jj2 [f{x, y) ■ g{x, y)f dx dy. 

Recall that Ojy is a value of ^ = {a,(3,d) = (a, /?, di, ^2) making a'j^iO) in 
(1.7) a minimum relative to the spectral density function f{x,y,9) given 
in (1.6). For simplicity we use 9 = {Oi, 02, 03,64) = {a, I3,di,d2). Define a 
symmetric matrix 

S = (aij) 

(3.1) f df-\x,y,eo) df-Hx,y,eo) 

with a., = j^^ / (x, y, 00) dx dy. 

A calculation shows that aij = fp ^ ^ sd%e '^°^ ' f{^^y^^o)dxdy and, more- 
over, using the spectral density function listed in (1.6), it is straightforward 
to verify 

(Til = 87r2/(l - ftp), 0-12 = 0, 

CT13 = -87r2(Log(l - ao))/ao, ctm = 0, 

(3.2) a22 = 8^V(l-/3o), fT23 = 0, 
cT24 = -87r2(Log(l-/3o))//3o, 

o"33 = (744 = 47r^/3 and 0-34 = 0. 

The interpretation given when oq = (/?o = 0) is that 0-13 (0-24) be replaced 
by its limiting value Svr^. Moreover, denote 

1 — 6 x^ 



and define 
K 



1 



87r2C(Qo)C(/?o)' 



7r2 



<5ii = — C(/?o), 5l2 = 0, 
D 

(3.3) Ji3 = ^°'^'""°^ C(/?o), <^i4 = 0, 
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h2 = -:rC(ao), 
b 



0, 5.4 = ^^^^^^C(ao), 



1 



^33 



l-«o 



cm, 



6m = and 5, 



1 



44 



C(ao). 



Then one can verify that S ^ = K{5ij). 

Proof of Theorem 1.2(a). Applying the mean value theorem, 
(3.4) 



for some On belonging to the line segment between and ^o- According 
to Theorem 1.1, 9]^ — > 9q almost surely and, thus, in probability; hence, 
-§e^Ni^N) = since becomes an interior point of the parameter space 
when is sufficiently large. It follows from Lemma 2.3, Theorem 1.1 and 
(3.1)-(3.2) that 

(3.5) An := ^^o'Ar(^Af) — > 5] almost surely as A^ — > oo. 

The Cramer-Wold device can be used to show that, as A^ — > oo, 



(3.6) 



N 



d_ 
do 



D 



AA(0,1287r6s). 



Indeed, let c' = (ci, C2, C3) G and define g{x, y, c) = c' 



d 



J df-H^,y,9o) . 



then 



1 



c' ■ N—aj^{9o) = N l_^g{x,y,c)lN{x,y)dxdy = —X'^B^Xn 



N 



where Bjy = T^ig)- According to Lemma 3.3, 



c'N 



d_ 

do 



^Nm-Ei^-^a%{9o)^ 



1 

N 



[X'j^BnXn - E{X'j^BnXn)] ^Af{0,S'') 



where 



6^ = 1287r6 / [fix, y, Oq) ■ g{x, y, c)f dx dy 

df-Hx,y,9o) df-\x,y,Oo) 



k,e=i 
1287r^ c'Sc. 



/2 



d9k 



dOt 



f{x,y,Oo)dxdy 
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Therefore, (3.6) holds. From this, employing (3.4) and (3.5), 

N{eM - On) + NA-}Ei —'2 



d 



de 



\d9 



D 



AA(0,1287r6s-i) 



as ^ oo. □ 



4. Proof of Theorem 1.2(b). According to (3.5), Aj^^ ^ S ^ almost 
surely and, thus, the proof of Theorem 1.2(b) involves determining the or- 
ders of NE{-§galf{eo)). Define aj^{9) as in (1.7) with Xjsi in /at replaced 
by /i = E{Xii). The argument given by Fox and Taqqu ([12], Lemma 8.1) 
verifies that NEl-^af^^Oo) — ^crj^iOo)] ^ as iV ^ co; it is shown that the 
orders of iV£;(^cr^*((9o)) and NE{-§galj{9o)) coincide. 

Assume that {Xst-s,t G Z} is a process described in Theorem 1.2 with 
spectral density function f(x,y,9) defined in (1.6) and having mean fj, and 
autocovariance function ^{k,£). Denote 



N 



kl=l 



and 



dxdy. 



1 12 f{x,y,B) 

It follows from the stationarity of the X-process that 



k,i,m,n=l 



X E{Xke - ij){Xynn - IJ) dx dy 



N 



iV2 



'y{k — m,£- 



n 



(4.1) 



k,£,m,n=l 



d 



X J^^ e^^^'^-^>+^'-^^y^-^f-\x,y,eo)dxdy 

E {N-mN-\i\h{k,e) 



\k\<N,\e\<N 



X J^^e'^'^+'y^^f-Hx,y,eo)dxdy. 
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Recall that 6q = (oq, (3o,dio,d2o) denotes the true parameters. The follow- 
ing notation is used: 

g{x,a,di) = |1 - e-^n-^'^ijl - 06"^^ "2, 
h{y,P,d2) = |1 - e-'y\-^''^\l - /3e~'y\-^ 



(4.2) 



flfci = J^e'''''g{x,ao,dio)dx, ai2 = J^e''^yh{y,Po,d2o)dy, 

e''^yh-^{y,Po,d2o)dy, Cki = J^e''''=-^{x,ao,dio)dx, 
f dg~^ 

i'iy) = a{y - ao, dio) "^^(^' "o, c?io) dx. 

It remains to determine the limit of the sequence {N {hog N)~^ E{^(jl^ {6q))} . 
For the sake of brevity, verification is presented here only for the di -component 
of 9. Employing the notation in (4.2), f~^{x, y, 6) = ^^g~^{x, a,di)- h~^{y, (3, ^2) 
and, thus, (4.1) becomes 



(4.3) E[±-aj^[e,))=Y. E (1 



fe|<7V ^ ' |^|<iV 



N 



£2- 



Observe that the second summation converges to Svr^ as ^ c« since the 
limit is precisely 47r^ times the Cesaro sum of the convolution of h and 
evaluated at zero. 



Lemma 4.1. Suppose that the assumptions listed in Theorem 2.1 are 
satisfied. Employing the notation of (4.2) gives the following: 

(i) V'(y) = -8diosin|Log^sin?/ + 0(y|Logy|) asy^Q+; 

(ii) Eii^aUOo)) = 0{^),E{f,al{eo)) = 0{^) and ^ - 



ddi 



E{^a%{9o))^-fTT^d^o,i = l,2, as 00. 



Proof. Verification of (i) involves applications of Lemmas A.l and A. 2 
in the Appendix, together with several technical arguments of approximation 
and expansion. The complete details of the proof are available on request 
from the authors. 

(ii). Proof of the third part for z = 1 is supplied here. According to (4.3), 
it remains to determine the limit of sequence ^^^^ ^ Vn '■= j-^^^ ^ S|fc|<Af(l — 

^)ofciCfci. Then Vn is the A^th Cesaro sum for ^ evaluated at zero and thus 



2i//2 ^ 
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has the integral representation Vat = Jjtp{y)K]\r{y) dy, where 

1 sin^Ny/l 

(4.4) KN{y) = { ' 

denotes the Fejer kernel (e.g., [7], page 71, or [26], page 88). Since ^ is an even 
function, ^Vn = J(^Q -i^/j^) 4){y)KN{y)dy + J^^^j^^^^i){y)KN{y) dy + 
I[c,7T]'^iy)^N{y)dy := Jm + Jn2 + ^Af3, where c e (0,1). Choosing c suf- 
ficiently small, according to Lemma A. 2, 'ip{y) = O(yLog^y) as y — > 0+ 
and, thus, | Jati] < KiN J^^Q^^^j^^yLog^ y dy = 0( ^°^^ ). Let HV'lli denote the 

Li-norm of V- Then \J^^\<i^\\^\\^=0{j^). 

It remains to estimate Jn2- Using the expansion for ip in (i) when c is 
sufficiently small and sufficiently large, 

Jn2= / ( -8(iiosin^Log^siny + 0(y|Logy|) )A'Ar(y)(iy. 

Jll/N,c) \ 2 J 

The integral of the second term is 0( ^°^^ ) and 
2 . sm^Ny/2 



f 2 . sin^Ay/2 , 
/ Log smy — jt^- dy 

J[i/N,c) smy/2 



Log^siny Log^ y\ . 2 , 
I ; 2 sm dy 

[i/N,c)\ smy/2 y J 2 

^ f ^2 sm'^Ny/2 , 
+ 2/ Log2y 

J[i/N,c) y 



f 0{y Log2 y) sin2 ^dy + 2 f Log^ y^J^^J^v/l ^y 
J[i/N,c) 2 Jii/N,c) y 

--■■h+h. 



The first integral /i is 0(1) and /2 = /[i/Ar,c) '^^-/[i/Ar.c) cos Nydy. 
Integrating by parts, it is shown that the second integral is 0(Log^ N) and, 
thus, l2 = l Log2 N + 0(Log2 A). Hence, 

f Log2 sin y^-^^^^m^ dy = h + h = 0(1) + ^ Log^ A + 0(Log2 A) 
Jli/N.c) smy/2 3 



^Log3A + 0(Log2A) 



and, thus. 



8dio /• ^ 2 • sin^Ay/2 /Log^A 
Jn2 = -7-^ / Log smy — ^ —dy + O' 



2ttNJ[i/n,c) smy/2 \ A 
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SvrTV \ N 



Combining the above results, 
Vn = 2( Jati + Jn2 + Jns) 



V N J \ 3ttN \ N J J \N 



SvriV V N 



and, thus, ^^^Vn as iV^ oo. It follows from (4.3) that 



T I l^'^N (^o) J ' — 5- 

Log^A^ \odi / Svr 3 

As mentioned earlier, NE{-^a'j^{6o) — ^o"^*(^o)) — ^ and, hence. 

Proof of Theorem 1.2(b). Recall from (3.3) and (3.5) that Aj/^ 
= K • {6ij) almost surely. According to Theorem 1.2(a) and Lemma 4.1, 

and NiLog^Nr^E{-§ga%{eo)) ^ where C = -^(0, 0, dio, 0(20). It fol- 
lows that 

N - P -1 

- — 3— -(6* AT -6*0) ^ -S 
Log A* 

In particular, A^(Log3 A^)~-^(Q!Ar — qq) ^ —K • J2^=i^ijCj = — -^^'^13^3 since 

6u = 0. Likewise, iV(Log3 N)~^0n - Po) ^ -i^524^4 since 623 = 0. Further, 

iV(Log3 iV)-i(di -dio) ^ -i^5336, iV(Log3 iV)-i('^^2 -f^2o) ^ -K6uU and, 
hence, Theorem 1.2(b) follows. □ 

APPENDIX 

The reader is referred to Boissy [6] for detailed proofs of the following 
lemmas. 



Lemma A.l. Assume that < d < 1 is fixed. Then for y > sufficiently 
small, 
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(i) < d(cos l)i-'^cot(^) sin I < ( Jg!(,)/2 )"-(cos f )'^< d(cos ^^f-'x 
cot(^-^)sin| when X G [—IT, —TT + y]; 

(ii) < -dcot(^)sinf < (cosf)'^ - U:^l%, r < -2dcot(2^)smf 
w/ien X G [— vr + y, — y] ; 

(iii) < d2'^-icot(2^)sin| < ( 4°.!g/2 )^ - < ^(cosD'^-i x 
cot(^^^) sin I when x £ [2y, vr] . 

Lemma A. 2. Fix < d < 1. Then there exist positive constants Ci,C2 
and 6 such that, for <y < 6, Ciy(Logy)^ < | |2sin ^|~'^|2sin fl*^ x 
Log 1 2 sin 1 1 dx\ < C2y(Logy)^. 
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